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The long-standing model-independent annual modulation effect measured by the DAMA Collab- 
oration, which fulfills all the requirements of a dark matter annual modulation signature, and the 
new result by the CoGeNT experiment that shows a similar behavior are comparatively examined 
under the hypothesis of a dark matter candidate particle interacting with the detectors' nuclei by 
a coherent elastic process. The ensuing physical regions in the plane of the dark matter-particle 
mass versus the dark matter-particle nucleon cross-section are derived for various galactic halo 
models and by taking into account the impact of various experimental uncertainties. It is shown 
that the DAMA and the CoGeNT regions agree well between each other and are well fitted by a 
supersymmetric model with light neutralinos which satisfies all available experimental constraints, 
including the most recent results from CMS and ATLAS at the CERN Large Hadron Collider. 



I. INTRODUCTION 

An annual modulation effect, as expected from the rel- 
ative motion of the Earth with respect to the relic parti- 
cles responsible for the Dark Matter (DM) in the galactic 
halo [H , has been measured by the DAMA Collaboration 
since long time with an increasing exposure along 13 
years which, with the second generation DAM A/LIBRA 
apparatus, has reached the value of 1.17 ton x year and 
a confidence level of 8.9 a 

A very recent analysis of the data collected by the 
CoGeNT experiment over a period of 442 days with a 
very low energy threshold Germanium detector having a 
fiducial mass of 330 g has now led this Collaboration to 
present the indication of a yearly signal modulation at 
about 2.86 a 0. 

The various experimental features required for a de- 
tector to be sensitive to the expected annual modulation 
effect are not met by most of the other direct detection 
experiments running at present. However, it is intrigu- 
ing that two of them (CDMS and CRESST 0) found 
in their data some excesses of events over what would be 
expected by them from backgrounds. It is also noticeable 
that, at least within one of the most widely considered 
kind of DM particles, i.e. the one with an clastic co- 



herent interaction with the atomic nuclei of the detector 
material, the CDMS and CRESST excess events would 
fall into (or close to) the physical region singled out by 
the DAMA/LIBRA and CoGeNT annual-modulation re- 
sults. 

The XENON100 Collaboration Q and the CDMS Col- 
laboration (in re-analyses of their previous data [§[) 
claim upper bounds - under a single set of fixed assump- 
tions - as in conflict with the aforementioned results of 
the other experiments. However, problems related to the 
conclusions of Refs. @, S], as discussed in Refs. [l3| 
and in Ref. 
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| , respectively, and the existence of many 
uncertainties in the procedures applied in the data han- 
dling by those experiments, lead us to carry out here an 
analysis of the results of Refs. 0,0] , n °t conditioned by 
the results reported in Refs. 0,13 - 

Though the model-independent annual modulation 
measured in the experiments of Refs. 0, § can be ac- 
counted for by a variety of interaction mechanisms of relic 
particles with the detectors materials [lij ]. we limit our 
analysis here to the case where the signal is caused by 
nuclear recoils induced by elastic coherent interactions 
with the DM particles. For simplicity as in a commonly 
used nomenclature, in the following we will call a generic 
particle with these features a WIMP (Weakly Interact- 
ing Massive Particle), although the term WIMP identi- 
fies a class of DM particles which can have well different 
phenomenologies, like e.g. a preferred interaction with 
electrons 131. 
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Thus, in Sect. [TTJ by using the results of Rcfs. 0,0] we 
first determine what are the physical regions pertaining 
to the DAMA and the CoGeNT annual modulation data 
in terms of the WIMP mass and of the WIMP-nucleon 
clastic cross-section at given confidence levels. In deriv- 
ing these regions we take into account the main origins 
of various experimental uncertainties, as well as differ- 
ent forms for the Distribution Function (DF) of DM relic 
particles in the galactic halo 14 1. 

Subsequently in Sect JIIII we show how the annual- 
modulation regions are well fitted by light neutrali- 
nos within the effective Minimal Supcrsymmctric exten- 
sion of the Standard Model (MSSM) at the electroweak 
(EW) scale introduced in Ref. (l5j . The relevance of 
light neutralinos in connection with the DAMA annual- 
modulation effect was first discussed in Ref. 16j; their 



phenomen olog y was then developed in the context of 
direct [17j, [l8[ and indirect [19J searches of DM parti- 
cles. The features of this specific realization of MSSM, 
dubbed Light Neutralino Model (LNM), are also con- 
fronted here with the most recent constraints on super- 
symmetry (SUSY) derived at the Tevatron and at the 
Large Hadron Collider (LHC). 

Particle physics models different from the LNM and 
potentially capable of generating light-WIMP particle 
candidates compatible with direct detection include su- 
persymmetric models which extend the MSSM by enlarg- 
ing the par ticle field content, like in the next-to-minimal 
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mirror- 



models [20j, sneutrino dark matter models 
dark matter models [22j], models with asymmetric dark 
matter (23|, isospin-violating models 24 1, singlet dark 
matter models 25j, specific realizations of grand uni- 
fication (26|, higgs-portal models 27 1, composite mod- 
els 



281 ] . specific two-higgs doublet models [29(; secluded 



31|. 



WIMPs [30[ • Additional recent analyses can be found in 
Ref. ' 



Conclusions of our analysis are drawn in Sect. [Vj 



II. REGIONS RELATED TO THE ANNUAL 
MODULATION EFFECT IN CASE OF WIMPS 

All experimental results discussed in the present Sec- 
tion are given in terms of plots in the plane m x - 



(nucleoli) 
scalar ' 



where ex. 



nucleoli) 
scalar 



is the WIMP-nucleon cross- 



section, £ = p x / po; po is the local total DM density and 
p x the local density of the DM candidate \. In the 
present section x denotes a generic WIMP candidate, 
main responsible for the annual modulation effect under 
discussion; it will specifically denote a neutralino in the 
Sections to follow. The factor £ leaves open the possibil- 
ity that the considered DM candidate does not provide 
the total amount of local DM density. 



A. Phase— space distribution functions of dark 
matter 

The quantity £&gcaXa* can ^ e derived from the exper- 
imental spectra, once a specific DF is selected to describe 
the phase-space distribution function of the DM particle 
in the galactic halo. The appropriate form for the DF is 
still the subject of extensive astrophysical investigation. 
It is also possible that DM direct detection might be af- 
fected by the presence of unvirialized components (see, 
for instance, Ref. \^ ) . Here we have taken a few samples 
of DFs, selected from the various realizations examined 
in Ref. 14j, specifically: i) the isothermal sphere (AO), 



ii) the Jaffe distribution ( A4) [33( , iii) a triaxial distribu- 
tion (D2) [53] (the notation adopted here follows those 
of Ref. [14j , to which we refer for further details) . Notice 
that one could also have DF with a non-isotropic veloc- 
ity dispersion (like distribution D2) and co-rotating or 
counter-rotating halos. Then, though the variety of DFs 
discussed in this paper already offer a significant sample 
of DFs, this selection is clearly not (and could not be) 
exhaustive of all possible situations. 

As for the main parameters characterizing the various 
DFs (the local total DM density po and the local rota- 
tional velocity vq) we will take into account their physical 
ranges as discussed in Ref. [13] ■ Thus, we will take as 
representative values of vq either one of the two extreme 
values or the central value of the physical range 170 km 
sec -1 < vo < 270 km sec _1 . For each representative 
value of i>o we take for po either its minimal po mm or its 
maximal value po max , i n the range compatible with the 
given value of Vq. As in Ref. [3] , po mm (po max ) is defined 
as the value to be associated to po when the visible mass 
provides its maximal (minimal) contribution to the total 
mass budget of the halo compatibly with observations. 
The numerical values for po min and po max depending on 
the DF and the values of vq will be taken from Table III 
of Ref. 



14j. The escape velocity will be set at v e 



=650 



km sec 



B. Annual-modulation regions in the considered 
model framework 



The about 9 a C.L. model independent positive re- 
sults of the DAMA/Nal and DAM A/LIBRA experiments 
and the recent positive hints by CoGeNT at 
2.86 a C.L. Q can be analyzed in many corollary model- 
dependent analyses. In all cases, many uncertainties 
in experimental parameters as well as in necessary as- 
sumptions on various related astrophysical, nuclear and 
particle-physics aspects must be taken into account. In 
the particular case of the WIMPs treated in this paper 
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many sources of uncertainties exist; some of them have 
been addressed e.g. in Refs. 35, 3^, 



. These affect all 
the results at various extent both in terms of exclusion 
plots and in terms of allowed regions/ volumes and thus 
comparisons with a fixed set of assumptions and param- 
eters values are intrinsically strongly uncertain. In the 
following we will point out the effect of just one exper- 
imental parameter, the quenching factor, whose precise 
determination is quite difficult for all kinds of used de- 
tectors. 

In fact, generally the direct measurements of quench- 
ing factors are performed with reference detectors, and 
- in some cases - with reference detectors with features 
quite different from the running conditions; in some other 
cases, these quenching factors are not even measured at 
all. Moreover, the real nature of these measurements 
and the used neutron beam/sources may not point out 
all the possible contributions or instead may cause uncer- 
tainties because e.g. of the presence of spurious effects 
due to interactions with dead materials as e.g. housing 
or cryogenic assembling, if any; therefore, they are intrin- 
sically more uncertain than generally derived. Thus, we 
specialize the present section to discuss the case of the 
values of the quenching factor of Na and I in the highly 
radiopure Nal(Tl) detectors of the DAMA experiments; 
analogous/similar discussions should be pursued for the 
other cases. 

As is widely known, the quenching factor is a specific 
property of the employed detector (s) and not a general 
quantity universal for a given material. For example, 
in liquid noble-gas detectors, it depends, among other 
things, on the level of trace contaminants which can 
vary in time and from one liquefaction process to an- 
other, on the cryogenic microscopic conditions, etc.; in 
bolometers it depends for instance on specific properties, 
trace contaminants, cryogenic conditions, etc. of each 
specific detector, while generally it is assumed exactly 
equal to unity. In scintillators, the quenching factor de- 
pends for example on the dopant concentration, on the 
growing method/procedures, on residual trace contami- 
nants, etc. and is expected to have energy dependence. 
Thus, all these aspects are already by themselves rele- 
vant sources of uncertainties when interpreting whatever 
result in terms of DM candidates inducing just recoils as 
those considered in the present paper. Similar ar gum ents 
have already been addressed e.g. in Refs. [i, 35, 36, 38j |. 
In the following, we will mention some arguments for 
the case of Nal(Tl), drawing the attention to the case 
of DAMA implications in the scenario considered in this 
paper. 

The values of the Na and I quenching factors used 
by DAMA in the corollary model-dependent calculations 
relative to candidates inducing just recoils had, as a first 



reference, the values measured in Ref. [40(. This mea- 
surement was performed with a small Nal(Tl) crystal ir- 
radiated by a 252 Cf source, by applying the same method 
previously employed in Ref. [4l|. Quenching factors 
equal to (0.4 ± 0.2) for Na and (0.05 ± 0.02) for I (in- 
tegrated over the 5 - 100 kcV and the 40 - 300 keV re- 
coil energy range, respectively) were obtained. Using the 



same parametrization as in Ref. 4l|, DAMA measured 
in Ref. [i^] quenching factors equal to 0.3 for Na and 0.09 
for I, integrated over the 6.5 — 97 keV and the 22 — 330 
keV recoil energy ranges, respectively. The associated 
errors derived from the data were quoted as one unity 
in the least significative digit. Then, considering also 
both the large variation available in the literature (see 
e^g. Table X of Ref. 



35 1 ) and the use of a test detector 



35j, a 20% associated error has been included. Never- 
theless, some recent considerations, as those reported in 
Ref. 44| about the energy dependence of quenching fac- 
tors for various recoiling ions in the same detector, have 
called our attention to the fact that the large uncertain- 
ties in the determination of Ref. 
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could be due, in a 
significant part, to uncertainties in the parametrization 
itself, which we also adopted. Another uncertainty could 
arise from the determination of integrated values, while 
an increase of the quenching factor values towards lower 
energies could be expected, as observed in some crystal 
detectors as for instance Csl. 

An additional argument on uncertainties on quench- 
ing factors in crystals, and specifically for Nal(Tl), is the 
presence and the amount of the well known channeling 
effect of low energy ions along the crystallographic axes 
and planes of Nal(Tl) crystals. Such an effect can have a 
significant impact in the corollary model dependent anal- 
yses, in addition to those uncertainties discussed above 
and later on, since a fraction of the recoil events would 
have a much larger quenching factor than that derived 
with neutron calibrations. Since the channeling effect 
cannot be generally put into evidence with neutron mea- 
surements, as discussed in details in Ref. 38J, only theo- 
retical modeling has been produced up to now. In partic- 
ular, the modeling of the channeling effect described by 
DAMA in Ref. [38[ is able to reproduce the recoil spec- 
trum measured at neutron beam by some other groups 
[39]] . For completeness, we mention alternative channel- 
ing models, as that of Ref. [42| where larger probabil- 
ities of the planar channeling are expected. Moreover, 
we mention the analytical calculation claiming that the 
channeling effect holds for recoils coming from outside a 
crystal and not from recoils produced inside it, due to the 
blocking effect 43|. Nevertheless, although some amount 
of blocking effect could be present, the precise description 
of the crystal lattice with dopant and trace contaminants 
is quite difficult and analytical calculations require some 
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simplifications which can affect the result. 

Recently, Ref. (3] pointed out the possibility that the 
quenching factors for nuclear recoils in scintillators can 
be described with a semi-empirical formula having only 
one free parameter: the Birks constant, fc^, which de- 
pends on the specific set-up. Applying this procedure 
to the DAMA detectors operating underground and fix- 
ing the ks parameter to the value able to reproduce the 
light response to alpha particles in these detectors, the 
expected Na and I quenching factors are established as a 
function of the energy with values ranging from 0.65 to 
0.55 and from 0.35 to 0.17 in the 2 - 100 kcV electron 
equivalent energy interval, for Na and I nuclear recoils, 
respectively; as evident, also an energy dependence is 
pointed out there. 

In the following analysis, we present some of the many 
possible model-dependent analyses of the DAMA results, 
including at least some of the present uncertainties. In 
particular, the uncertainties due to the description of the 
halo are accounted for some of the many possible halo 
models; we employ here the DFs mentioned in Sect III Al 



(nucleoli) 
scalar 



as a function of the 



Figs. [U H and E] show £cr 
Dark Matter particle mass m x for the AO, A4 and D2 
halo models 14[. In order to have a significative sample 
in terms of the physical ranges of the relevant astrophys- 
ical parameters we have chosen to display the annual- 
modulation regions for the two extreme values of the local 
rotational velocity vq, i.e.: vq = 170 km sec -1 (in the left 
panel of each figure) and v = 270 km sec -1 (in the right 
panel). In Fig. [TJ where the case of the AO distribution 
function is shown, we have set the local total DM density 
po to be equal to its minimal value, po = po mm , compat- 
ibly with the value of vq, then po = 0.18 GeV cm~ 3 for 
vq = 170 km sec -1 (in the left panel) and po = 0.45 GeV 
cm for Vq = 270 km sec -1 (in the right panel). In Fig. 
[21 where we display the case of the A4 distribution func- 
tion, we set again p = po min , then p = 0.26 GeV enr 3 
for vq = 170 km sec -1 (in the left panel) and po = 0.66 
GeV cm -3 for vo = 270 km sec -1 (in the right panel). 
Fig. [3] shows the case for the D2 distribution function 
for a value of po equal to its maximal value, po = Po max ; 
thus, po = 0.50 GeV cm" 3 for vq = 170 km sec -1 (in 
the left panel) and po = 1.27 GeV cm -3 for vq = 270 
km sec -1 (in the right panel). The values for po mm and 
po max employed here are taken from Table III of Ref. 
14 |. A further example of annual-modulation regions, 



corresponding to the standard DF, the cored-isothermal 
sphere AO with p = 0.34 GeV cm -3 and v = 220 km 
sec -1 , will be given in Fig. [6] in Sect. IIV1 where we 
compare the LNM with the experimental results. In all 
figures the escape velocity has been maintained at the 
fixed value: 650 km/s. Of course, the present existing 
uncertainties affecting the knowledge of the escape ve- 



locity - as well as other uncertainties not included here 
- would significantly modify/extend the allowed regions. 

The three (colored) hatched regions in Figs. [TJ [2] and 
[3] denote the DAMA annual modulation regions, un- 
der the hypothesis that the effect is due to a WIMP 
with a coherent interaction with nuclei and in 3 differ- 
ent instances: i) without including the channeling ef- 
fect ((green) vertically-hatched region), ii) by includ- 
ing the channeling effect according to Ref. [38( ((blue) 
horizontally-hatched region), and iii) without the chan- 
neling effect but using an energy-dependent Na and I 
quenching factors as established by the procedure given 
in Ref. [44| ((red) cross-hatched region). It is worth not- 
ing that, depending on the possible amount of blocking 
effect in Nal(Tl) with respect to the modeling used in 
Ref. [i^], the channeled (blue) region will span the do- 
main between the present channeled region and the un- 
channeled one. Moreover, the availability of quenching 
factor values not integrated over a large energy interval 
can also play a relevant role. 

All these DAMA regions have been investigated here in 
some specific cases by adopting a procedure that allows 
to put into evidence - to some extent - the uncertainties 
on the quenching factors and on the nuclear form fac- 
tors: by considering the mean values of the parameters 
of the used nuclear form factors and of the quenching 
factors of Ref. [i^] (case A); by varying the mean val- 
ues of those quenching factors up to +2 times the errors 
quoted there and the nuclear radius, r n , and the nuclear 
surface thickness parameter, s, in the SI form factor from 
their central values down to —20% (case B); by fixing the 
Iodine nucleus parameters at the values of case B, while 
for the Sodium nucleus one considers the quenching fac- 
tor at the lowest value measured in the literature and the 
nuclear radius, r„, and the nuclear surface thickness pa- 
rameter, s, in the SI form factor from their central values 
up to +20% (case C). 

The DAMA regions have been obtained by superposi- 
tion of the three regions corresponding to the cases A, 
B, and C. These regions represent - as in some previous 
DAMA publications - the domain where the likelihood- 
function values differ more than 7.5 a from the null hy- 
pothesis (absence of modulation). This choice allows 
both a direct superposition of the obtained results for 
both Na and I target nuclei (which case by case can have 
different levels of the corresponding minimum value of 
the likelihood function) and a very high C.L. require- 
ment. 

In the same figures Figs lll2l3l the allowed regions by 
the CoGeNT experiment [4| (under the same adopted 
framework) are reported, assuming for simplicity for the 
Ge a fixed value of 0.2 for the quenching factor and a 
Helm form factor with fixed parameters. In particular, 
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m x (GeV) m x (GeV) 

FIG. 1: ^Pscalai as a function of the mass m x of a generic DM particle which interacts with nuclei by an elastic coherent 
scattering. The halo DF is taken to be given by the isothermal sphere ((AO) in the notations of Subsect lll Al and Ref . ) . The 
parameters are: i) in the left panel, wo = 170 km sec -1 , po = 0.18 GeV cm -3 ii) in the right panel, wo = 270 km sec -1 , po = 0.45 
GeV cm -3 (see text for further details). The three (colored) hatched regions denote the DAMA annual modulation regions, 
under the hypothesis that the effect is due to a WIMP with a coherent interaction with nuclei and in 3 different instances: 
i) without including the channeling effect ((green) vertically-hatched region), ii) by including the channeling effect according 
to Ref. [3^] ((blue) horizontally-hatched region), and iii) without the channeling effect but using the energy-dependent Na 
and I quenching factors as established by the procedure given in Ref. [3] ((red) cross-hatched region). They represent the 
domain where the likelihood-function values differ more than 7.5 a from the null hypothesis (absence of modulation). The 
(non-hatched) region denoted by a (black) solid contour is the allowed region by the CoGeNT experiment when considering 
the modulation result given in Ref. and the assumptions given in the text for the quenching factor and the form factor. 
This region is meant to include configurations whose likelihood-function values differ more than 1.64 a from the null hypothesis 
(absence of modulation). This corresponds roughly to 90% CL far from zero signal. In fact due to the presently more modest 
C.L. (about 2.9 a) of this result with respect to the 9 a C.L. of the DAMA/Nal and DAMA/LIBRA evidence for Dark Matter 
particles in the galactic halo, no region is found if the stringent 7.5 a from absence of modulation is required as for DAMA. It 
is worth noting that, depending on other possible uncertainties not included here, the channeled (blue) region could span the 
domain between the present channeled region and the unchanneled one. 



the CoGeNT regions have been obtained by fitting the 
measured modulation amplitudes with the WIMP expec- 
tation (S m ) and using the 0.45-3.15 keV energy region (R 
in the following) of the energy spectrum as a constraint. 
The x 2 function is: 

x = }^ — -2 + 1^ — -2 e(s* , fc -r fe ), 

fc=l,2 R 

(1) 

where A k and &A,k are the modulation amplitudes and 
their errors in the two considered energy bins; r k and 
(Tfc are the rates and their errors in the k energy bin. 
The Heavisidc function occurs in the second term 
to account for the constraint of the rate in those en- 
ergy bins (R). In particular, we derived from Ref. 
the following modulation amplitudes: 4(0.5 — 0.9) 
keV = (0.91 ± 0.61) cpd/kg/kcV; 4(0.5 - 3.0) keV = 



(0.45 ± 0.18) cpd/kg/keV. Thus, we consider in eq. Q] 
4 fc=1 = .4(0.5-0.9) keV and we infer A h=2 = A(0.9-3.0) 
keV = (0.36 ±0.18) cpd/kg/keV. The values of the mod- 
ulation amplitudes have been obtained here under the 
assumption that the period and the phase of the mod- 
ulation are fixed at their nomimal values of 1 year and 
June 2 nd . If one allows the phase and the period to be 
free parameters, the ensuing modulation amplitudes oc- 
cur to be larger, but still compatible within the quoted 
errors. 

The (non-hatched) regions denoted by (black) solid 
contours in Figs. [TJ [2] and [3] denote the allowed regions 
by the CoGeNT experiment; such regions contain con- 
figurations whose likelihood-function values differ more 
than 1.64 a from the null hypothesis (absence of mod- 
ulation). This corresponds roughly to 90% CL far from 
zero signal. In fact due to the presently more modest 
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FIG. 2: As in Fig. [T] except that here the halo DF is taken to be given by the Jaffe distribution ((A4) in the notations of 
Subsect III Al and Ref.[3l). The parameters are: i) in the left panel, vo = 170 km sec -1 , po = 0.26 GeV cm -3 ii) in the right 
panel, vo = 270 km sec" 1 , po = 0.66 GeV cm -3 




FIG. 3: As in Fig. [T]except that here the halo DF is taken to be given by a triaxial distribution [3J] ((D2) in the notations of 
Subsect III Al and Ref.[3]). The parameters are: i) in the left panel, no = 170 km sec -1 , po = 0.50 GeV cm -3 ii) in the right 
panel, vo = 270 km sec -1 , po = 1.27 GeV cm -3 



C.L. (about 2.9 a) of the CoGeNT result with respect to 
the 9 g C.L. of the DAMA/Nal and DAM A/LIBRA evi- 
dence for DM particles in the galactic halo, obviously no 
region is found if the stringent 7.5 a from absence of mod- 
ulation is required as for the DAM A cases; thus, it will 
be very interesting to see future CoGeNT data releases 



with increased significance. Anyhow, all the examples 
given here, as well as the proper inclusion of possible 
uncertainties in the assumptions adopted for CoGeNT, 
and additional accounting of other uncertainties, offer a 
substantial agreement between the two experiments (as 
well as with some preliminary possible positive hint by 
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CRESST discussed at Conferences so far Q, which is not 
addressed here) towards a low mass candidate. 

From Figs. [TJ [2] and [3] we see that in all instances the 
DAM A and the CoGeNT regions agree quite well, over 
ranges of C^caiar* an d m x somewhat wider as com- 
pared to those derived for instance in Refs. 0, |45|. The 
gross features in the comparative positions of the vari- 
ous regions in our Figs. [TJ [5] and [3J are easily understood 
in terms of the specific values of the DF parameters em- 
ployed. Further statistics in the CoGeNT experiment will 
be useful in pinning down more precisely the common do- 
mains for the two annual-modulation experiments. 



III. THE LIGHT NEUTRALINO MODEL 

Now we discuss how the results reported in the previ- 
ous Section arc well fitted by the light neutralinos which 
arise within the model introduced in Ref. 15j and de- 
veloped in the papers of Ref. Lately, this model, 



denoted as Light Neutralino Model (LNM), was updated 



in Refs. [18j, |46[ to take into account recent constraints 
on supersymmetric models derived at accelerators and 
B-factories. 



Since no gaugino-mass unification at a Grand Unified 
scale is assumed in our LNM (at variance with one of 
the major assumptions in mSUGRA), in this model the 
neutralino mass is not bounded by the lower limit m x > 
50 GeV that is commonly derived in mSUGRA schemes 
from the LEP lower bound on the chargino mass (of 
about 100 GeV). In Refs. [Il-QII it is shown that, if R- 
parity is conserved, a light neutralino (i. e. a neutralino 
with m x < 50 GeV) is a very interesting candidate for 
cold dark matter (CDM), due to its relic abundance and 
its relevance in the interpretation of current experiments 
of search for relic particles; it is also shown there that a 
lower bound m x > 7-8 GeV is obtained from the cosmo- 
logical upper limit on CDM. The compatibility of these 
results with all experimental searches for direct or indi- 
rect evidence of SUSY (prior to the first physics results of 
LHC) and with other precision data that set constraints 
on possible effects due to supcrsymmetry is discussed in 
detail in Ref. 18j . The viability of very light neutralinos 
in terms of various constraints from collider data, preci- 
sion observables and rare meson decays is also considered 
in Ref. 59|. Perspectives for investigation of these neu- 
tralinos at LHC are analyzed in Ref. and prospects 
for a very accurate mass measurement at ILC in Ref. 
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Main features of the LNM 



The LNM is an effective MSSM scheme at the elec- 
troweak scale with the following independent parame- 
ters: Mi, M2, M3, fi, tan/3, tua, tnq, ttij and A. Notations 
are as follows: Mi, M 2 and M 3 are the U(l), SU(2) 
and SU(3) gaugino masses (these parameters are taken 
here to be positive), /1 is the Higgs mixing mass param- 
eter, tan/3 the ratio of the two Higgs v.e.v.'s, tua the 
mass of the CP-odd neutral Higgs boson, rriq is a squark 
soft-mass common to the all families, mj is a slepton 
soft-mass common to all sleptons, and A is a common 
dimensionless trilinear parameter for the third family, 
A~ h = Ai = Arriq and Af = Amj (the trilinear parameters 
for the other families being set equal to zero). We recall 
that in Ref. [46| the possibility of a splitting between the 
squark soft-mass common to the first two families and 
that of the third family was considered. This allows to 
reduce the fine-tuning in the parameters that can be in- 
duced by the interplay between the constraint from the 
b —> S7 decay and those from SUSY searches at the LHC. 

The linear superposition of bino B, wino and of 

the two Higgsino states H^, H% which defines the neu- 
tralino state of lowest mass m x will be written here as: 



In the present section we essentially recall the main 
properties of the light neutralinos within the LNM, as 



derived in Refs. 15-181, and which arc relevant for the 



X = aiB + a 2 W {3) + a 3 H^ + a A H% 



(2) 



discussion of the experimental results of Refs. 



In the regime of light neutralinos the lower limit on 
the mass m x , obtained from the requirement that its relic 
abundance does not exceed the observed upper bound for 
cold dark matter (CDM), i.e. tt x h 2 < (flc DM^^max, 
can be expressed analytically in terms of the relevant 
SUSY parameters. In this concern, it is convenient to 
distinguish between two scenarios. The first one is de- 
noted as Scenario A, and its main features are: i) niA is 
light, 90 GeV < m A < (200 - 300) GeV (90 GeV being 
the lower bound from LEP searches); ii) tan/3 is large: 
tan j3 = 20-45, iii) the B—H° mixing needs to be sizeable, 
which in turn implies small values of yu: ~ (100 — 200) 
GeV. In this scenario the dominant contribution to the 
annihilation cross-section of a pair of neutralinos, <7 ann 
(which establishes the size of the neutralino relic abun- 
dance) is provided by the j4-exchange in the s channel 
of the annihilation process \ + X ~ * b + b, thus the lower 
bound on m x is given by: 



8 



m . 



[1 - m 2 b /ml 



(2m, 



) 2 /m\ 



> 7.4 GcV 



35 
^tan/3 

I 



0.12 

2 2 



0.12 



1 



(3) 



where mj, is the mass of the b quark. 

When niA > (200 — 300) GeV, the cosmological upper 
bound on the neutralino relic abundance can be satisfied 
by a pair annihilation process which proceeds through 
an efficient stau-exchange contribution (in the t, u chan- 
nels). This requires that: (i) the stau mass rrif is suffi- 
ciently light, rrif ~ 90 GeV (notice that the current ex- 
perimental limit is rrif ~ 87 GeV) and (ii) \ is a very 
pure Bino (i.e. (1 - a\) ~ 0(1O~ 2 )). Thus, one is 



lead to a Scenario £>, identified by the following sector 
of the supersymmetric parameter space: Mi ~ 25 GeV, 
\fi\ > 500 GcV, tan/3 < 10; m r > (100 - 200) GcV, 
-2.5 < A < +2.5. As derived in Ref. [H-QJ], in this 
scenario the cosmological bound fl x h 2 <J£lcDMh 2 ) m ax 
provides the lower bound m x > 22 GeV [47|, whose scal- 
ing law in terms of the stau mass and (ficDM^Jmax is 
approximately given by: 



J 



m 



: [l-ml/m 2 x } 1/4 > 22 GeV ( 



90 GeV 
I 



0.12 



r 



(4) 



In Scenario A the neutralino-nuclcon cross-section h, whose mass m/, has a numerical value very close to 



(nucleoli) 
scalar 



is dominated by the interaction process due to 
the exchange of the lighter CP-even neutral Higgs boson 



(nucleoli) . 

niA', then er , is expressible as: 



J 



^° n) - 5.3 x 10- 



■41 2 

cm 



9 

0.13 



tan^\ 2 /90 GcV\ 4 / g d 



35 



r 



V m h 



290 McV 



(5) 



where gd is the dominant coupling in the interaction of 
the Higgs boson h with the d-type quarks, 



9d = [m d (N\dd\N) + m s (N\ss\N) + m b (N\bb\N}}, (6) 

and (N\dd\N) denotes the scalar density of a generic 
quark q inside the nucleon. In this expression we have 
used as reference value for g d the v alue gd,ref 

= 290 MeV 



employed in our previous papers [17} . We recall that 
this quantity is affected by large uncertainties [48| with 

(ffd.max/ffd.ref) 2 = 3.0 and (gd,min/ 9d,rei f = 0.12, a fact 

that directly transforms in the same amount of uncer- 
tainty on the coherent scattering cross section. 

Since, as mentioned in Sect. UH we wish to con- 



sider also situations where relic ncutralinos only provide 
a fraction of the CDM abundance, the relevant quan- 
tity we will compare with the experimental results is 

i (nucleon) -i , , n j. (nucleon) mi r 

not simply a y scalar ', but rather ^< alar ; . The fac- 
tor £ = p x / po is calculated here according to the rescal- 



min{l,n x /i 2 /(f2 C DMfr 2 )min} (49j, where 



ing recipe £ 

(^c.DM^ 2 )min is the minimal value to be assigned to the 
relic abundance of CDM. 

It is remarkable that for neutralino configurations, 
whose relic abundance stays in the cosmological range 
for CDM (i.e. (VLcDAih 2 ) mi „ < fl x h 2 < (flcDAih 2 )^ 
with (£lcDAfh 2 ) m in = 0.098 and (0 C DJi/'i 2 )max = 0.12) 
and pass all particle-physics constraints, the elastic 
neutralino-nucleon cross-section can be cast as 



J 



(nucleon) 
^scalar 



(2.7-3.4) x 10 



41 cm 2 



( 9d V 

V 290 McV/ 



[1 - (2m x ) 2 /m 



2 12 



(m x /(10 GeV) 2 [1 



rn^/m 2 ] 1 / 2 



(7) 



Notice that this formula provides an evaluation of 
CT scaiai COn ^ simply m terms of the neutralino mass m x . 
Specific SUSY parameters such as tan/3 or do not ap- 
pear explicitly, since these parameters have been reab- 
sorbed by the introduction of the relic abundance (this is 
because here the annihilation amplitude is related to the 
clastic-scattering amplitude by crossing symmetry) . The 
remaining dependence on the mass of the interaction me- 
diator rriA — mh is only marginal, due to the small values 
of m x considered here. 

Eq. {7} is of particular interest in establishing the range 
of values for cr^^ OIL ^ m terms of the neutralino mass. 
The numerical range in front of Eq.Q follows from the 
requirement that relic neutralinos have an abundance in 
the cosmological range for CDM. The crucial factor of 
uncertainties in Cg""^ 0011 ' is related to QCD properties 
through the coupling g^. It is however worth recalling 
that the range of the neutralino mass depends on the 
lower bound on m x which is explicitly given in terms of 
the SUSY parameters in Eq.®. These properties will 
also show up later in the figures displaying the scatter 

-i , r {nucleoli) 

plots for ^ calar ; . 



B. Constraints on SUSY parameters from early 
searches at the LHC 
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FIG. 4: Upper bounds in the mA - tan /3 plane, derived in 
Ref. from searches of the neutral Higgs boson decaying 
into a tau pair at LHC [5^|. The disallowed domain is the 
(yellow) shaded region. The solid bold lines labeled by num- 
bers denote the cosmological bound Q. x h 2 < (flcDMh 2 )^^ 
for a neutralino whose mass is given by the corresponding 
number (in units of GeV), as obtained by Eq. © with 
(QcDMh 2 )mnx = 0.12. For any given neutralino mass, the 
allowed region is above the corresponding line. 



In Ref. [46| the possible impact of some early analyses 
by the CMS and ATLAS Collaborations at the LHC on 
the LNM was investigated. The data considered there 
consisted in the results of searches for supersymmetry 
in proton-proton collisions at a center-of-mass en ergy 
of 7 TeV with an integrated luminosity of 35 pb~ [50(, 
i.e. the results of the CMS Collaboration for events with 



neutralino mass. It was proved there that the data of 
5ol . 51 1 do not imply a modification of the lower 



jets and missing transverse energy [50| . and those of the 
ATLAS Collaboration by studying final states containing 
jets, missing transverse energy, either with a isolated lep - 
ton (electron or muon) [sil j or without final leptons [H^j. 
Both signatures would be significant of processes due to 
the production in pairs of squarks and gluinos, subse- 
quently decaying into quarks, gluons, other standard- 
model (SM) particles and a neutralino (interpreted as 
the lightest supersymmetric particle (LSP)) in a R-parity 
conserving SUSY theory. As reported in Refs. [5fJ, |51| 
the data appeared to be consistent with the expected 
SM backgrounds; thus lower bounds were derived on the 
squark and gluino masses which are sizeably higher than 
the previous limits established by the experiments DO 
and CDF at the Tevatron. 
These data were employed in Ref. 



46] to determine 



the relevant lower bounds on the squark masses and the 
gluino mass M3 within the LNM and their ensuing pos- 
sible impact on the value of the lower bound on the 



Refs. 

bound m x > 7-8 GeV for the LNM, when the common 
squark mass for the first two families rriq 12 and the one 
for the third family mj are independent parameters with 
m,q 12 > mj, or, in case of a full degeneracy of the squark 
masses over the 3 families (as considered in the present 
paper), when M3 > (1.5 - 2) TeV. Otherwise, in the case 
of a full squark-mass degeneracy (rnq 12 = mi = m^) the 
lower bound on m x varies as a function of the gluino 
mass M 3 , from the value of 7-8 GeV for M 3 > 2 TeV 
to about 12 GeV for M 3 ~ 600 GeV (see Fig. 5 of Ref. 

for details). In particular, the gluino mass enters in 
the calculation of observables for the relic neutralino only 
at the loop level (through radiative corrections of Higgs 
couplings), so within the LNM M3 is very weakly cor- 
related to the other parameters. In order to reduce the 
number of parameters, in the present analysis we choose 
to decouple the gluino mass assuming M3 =2 TeV. In 
this case LHC data imply the lower bound m^ > 450 
within the LNM [3]. In the following we will impose 
this constraint in our numerical analysis. 

Now, we proceed to a discussion of the new results pre- 
sented by the CMS Collaboration on a search for neutral 
SUSY Higgs bosons decaying in tau pairs at a center- 
of-mass energy of 7 TeV with an integrated luminosity 
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FIG. 5: Scatter plot for f^^— 6 as a function of the neutralino mass for <?d,ref = 290 MeV. The left panel displays by (red) 
crosses SUSY configurations with a neutralino relic abundance which matches the WMAP cold dark matter amount (0.098 
< Q, x h 2 < 0.122) whereas the right panel displays by (blue) dots the configurations where the neutralino is subdominant 
(fl x h 2 < 0.098). The (light-blue) flag-like region denotes the extension of the scatter plot upwards and downwards, when the 
hadronic uncertainties are included. 



of 36 pb -i |55|. Since no excess is observed in the tau- 
pair invariant -mass spectrum, upper limits on the Higgs- 
boson production cross section times the branching ratio 
to tau pairs are placed. These limits are then converted 
into upper bounds for the SUSY parameter tan (3 as a 
function of the pseudoscalar Higgs-boson mass ttla in a 
particular MSSM benchmark. The ensuing disallowed re- 
gion in the plane ttla — tan /3 turns out to be considerably 
larger than the one previously derived at the Tevatron 
(see for instance Ref. [56| ) . 

However, in Ref. [H3] it has been shown that, when all 
the theoretical uncertainties involved in the derivation 
of the previous bounds on the Higgs-boson production 
cross section times the branching ratio to tau pairs are 
appropriately taken into account, the limits on the SUSY 
parameters reported in Rcfs. 
relaxed. 

We display in Fig. [4] the region disallowed in the plane 
(tan (3 — ijia) from the results of Rcfs. 
the analysis of Ref. 



55l . |56| are significantly 



55j , as derived in 



57j . In this figure we also show the 



lines corresponding to fixed values of the neutralino mass 
in the LNM. Thus we see that the CMS upper bounds 
of Ref. [55| do not modify the value of the neutralino- 
mass lower bound m x > 7-8 GeV, previously derived in 
Refs. 17, 18[. This result also follows directly from the 



analytic expression of Eq. (|3|) for the lower limit on m x . 

The predictions of the LNM for the cross-section 
^ were already anticipated in the analytic ex- 
Now we give the numerical 
which will finally be 



scat 



pressions of Eqs.([I)])-([7]). 
values for the quantity £ CT scaiar °^ 



confronted with the experimental results. Fig. [5] pro- 
vides the scatter plots of this quantity for the neutralino 
configurations which pass all the constraints previously 
discussed in this Section. In particular, in our scan of 
the LNM the following ranges of the parameters are 
adopted: 10 < tan/3 < 50, 105 GeV < n < 1000 GeV, 
5 GeV < Mi < 50 GeV, 100 GeV < M 2 < 2500 GeV, 
450 GeV < mg < 3000 GeV, 115 GeV < mj < 3000 GeV, 
90 GeV < m A < 1000 GeV, -3 < A < 3. 

The left panel refers to SUSY configurations with a 
neutralino relic abundance which matches the WMAP 
cold dark matter amount (0.098 < Q x h 2 < 0.122) 
whereas the right panel displays by (blue) dots the con- 
figurations where the neutralino is subdominant (fl x h 2 < 
0.098). In both panels, the flag-like region denotes the 
extension of the scatter plots upwards and downwards, 
when the hadronic uncertainties are included. 



IV. COMPARISON WITH EXPERIMENTAL 
RESULTS 

The predictions for C^caiar 00 "'' f° r light neutralinos 
within the LNM, as depicted in Fig. fall clearly in 
the region of interest of the present annual-modulation 
results as reproduced in Figs. I1I2I3I and in Fig. [7] (to 
follow). For a more specific comparison among experi- 
ments and theory we employ here, as a reference DF, the 
standard isothermal sphere with parameters: po = 0.34 
GeV cm~ 3 , vo = 220 km sec -1 , v esc = 650 km sec -1 . 
This choice is not meant to attribute to this particular 




m (GeV) 

FIG. 6: tj' 11 " 01 ™' as a function of the neutralino mass. The 

~ SCcllclr 

experimental annual-modulation regions are obtained as ex- 
plained in the caption of Fig.l, except that here the used DF is 
an isothermal sphere with the following values for the parame- 
ters: po — 0.34 GeV cm -3 , vo = 220 km sec -1 , v esc = 650 km 
sec -1 . The theoretical scatter plot displays the whole sample 
of neutralino configurations: (red) crosses denote SUSY con- 
figurations with a neutralino relic abundance which matches 
the WMAP cold dark matter amount (0.098 < Q x /i 2 < 0.122) 
while (blue) dots denote the configurations where the neu- 
tralino is subdominant (Q x h 2 < 0.098) (these two sets of con- 
figurations were shown separately in Fig.5). The scatter plot 
has been evaluated for gd,rd = 290 MeV. The (light-blue) 
flag-like region denotes the extension of the scatter plot up- 
wards and downwards, when the hadronic uncertainties are 
included (see text). 



DF a privileged role over other DFs, but is done sim- 
ply for convenience, to conform to the most commonly 
employed form of the DF. The experiment-theory com- 
parison is therefore displayed in Fig. [51 The features 
of this figure confirm that the conclusions drawn in Rcf. 
[l8j are even reinforced, when, as done in the present pa- 
per, the new annual-modulation results by CoGeNT are 
included; specifically: i) the light neutralino population 
agrees with the DAM A/LIBRA annual modulation data 
over a wide range of light neutralinos: 7-8 GeV < m x < 
50 GeV, ii) this population is also in agreement with the 
data of CoGeNT in a range of the neutralino mass some- 
what restricted to the lower masses: 7-8 GeV < m x < 
(15-20) GeV. 

It is worth recalling that also the data of CDMS 
and CRESST should their reported excesses be sig- 
nificant of real DM signals, would fall in a domain of the 
^scalar ^ ~ m x P iane overlapping with the DAMA and 
CoGeNT regions. 
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FIG. 7: Regions in the ^CTg^^° vs m x plane allowed by 
DAMA experiments in the three considered instances for the 
Na and Iquenching factors, including all the DFs considered 
in Ref. [3 and the same uncertainties as in Refs. [35LI3H] for 
a WIMP with a pure SI coupling. The hatchings (and colours) 
of the allowed regions are the same as those in Fig. [TJ These 
regions represent the domain where the likelihood-function 
values differ more than 7.5 a from the null hypothesis (ab- 
sence of modulation). It is worth noting that, depending on 
other possible uncertainties not included here, the channeled 
(blue) horizontally-hatched region could span the domain be- 
tween the present channeled region and the unchanneled one. 
The allowed region obtained for the CoGeNT experiment, in- 
cluding the same astrophysical models as in Ref. [35I . [3r3 [ and 
assuming for simplicity a fixed value for the Ge quenching 
factor and a Helm form factor with fixed parameters, is also 
reported and denoted by a (black) thick solid line. This region 
is meant to include configurations whose likelihood-function 
values differ more than 1.64 a from the null hypothesis (ab- 
sence of modulation). This corresponds roughly to 90% CL 
far from zero signal. See text. 

V. CONCLUSIONS 

The long-standing model-independent annual modu- 
lation signal measured by the DAMA Collaboration for a 
total exposure of 1.17 ton x year and a confidence level of 
8.9 a with a Nal(Tl) detector Q has been comparatively 
examined with the new results by the CoGeNT experi- 
ment [3| which shows a similar behavior with a statisti- 
cal significance of about 2.86 a. The annual modulation 
measured in these two experiments is an effect expected 
because of the relative motion of the Earth with respect 
to the relic particles responsible for the dark matter in 
the galactic halo 0. The underlying physical process 
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can be due to a variety of interaction mechanisms of relic 
particles with the detector materials [l2|- Here we have 
limited our analysis to the case where the signal would 
be caused by nuclear recoils induced by elastic coherent 
interactions of the target nuclei with the DM particles. 

The ensuing physical regions in the plane of the DM- 
particle mass versus the DM-particle - nucleon cross- 
section have been derived for a variety of DM distribution 
functions in the galactic halo and by taking into account 
the impact of various experimental uncertainties. 

The phase-space distribution of DM particles in the 
halo is still subject of extensive astrophysical investiga- 
tions, with the possible presence of unvirialized compo- 
nents (sec, for instance, Ref. [l^l). Here we have selected 
a few samples of DFs, selected among those discussed in 
Ref. [3] , from the isothermal sphere to the Jaffe DF 
to a triaxial one 



34|. 



We have examined in details to what extent the major 
experimental uncertainties, most notably those related to 
the quenching factors and the channeling effect, affect the 
derivation of the annual-modulation physical regions. It 
is shown that the DAMA and the CoGcNT regions agree 
well between each other independently of the specific an- 
alytic form of the DFs considered here, considering also 
that some existing uncertainties have not been taken into 
account for the CoGeNT allowed regions. For complete- 
ness, Fig. [7] shows the DAMA allowed regions in the three 
considered instances for the Na and I quenchin g fa ctors 
when including all the DFs considered in Ref. 14[ and 
the same uncertainties as in Ref. 3^, 36|. The allowed 



region obtained for the CoGeNT experiment, including 



the same astrophysical models as in Ref. [35|, |36[ and as- 
suming for simplicity a fixed value for the Ge quenching 
factor and a Helm form factor with fixed parameters, is 
also reported in Fig. [7J (solid line); it fully overlaps the 
DAMA allowed regions. The inclusion of other uncer- 
tainties on parameters and models (see for example Refs. 



35 



36j ) would further enlarge these regions. 
In this paper, we have finally discussed a specific 
particle-physics realization, the Light Neutralino Model, 
where neutralinos with masses in the tens of GeV range 
naturally arise. This supersymmetric model, which was 
already shown [l~8| to be successful in fitting the DAMA 
annual modulation results Q as well as the (unmodu- 
lated) CoGeNT (ill, the CDMS and the CRESST @ 
excesses, is shown here to agree quite well also with the 
most recent CoGeNT annual-modulation data 0. No- 
tice that the LNM discussed here satisfies all available 
experimental particle-physics constraints, including the 
most recent results from CMS and ATLAS at the CERN 
Large Hadron Collider. Confirmation of the validity of 
the SUSY model discussed in the present paper rests on 
the possibility of a positive evidence of light neutralinos 



in further running of LHC [6 Of . 

Note Added. In this Note we comment on two 
preprints that appeared after the submission of the 
present paper. 

The viability of an MSSM to obtain a neutralino- 
nucleon elastic cross section with a size relevant for the 
DAMA and CoGeNT annual modulation data, for light 
neutralino masses, is questioned in the preprint of Ref. 
62 1 . We note that in the MSSM scheme employed in Ref. 
62j the squark masses are all set at 1 TeV. From the 



properties discussed in detail in Refs. [18|, |46j it is clear 
that, taking all the squark masses at this value generates 
tension between the b — > s+7 and the B s — > con- 
straints, and thus precludes low values of the Higgs-boson 
masses (i.e close the their LEP lower bounds). This, in 
turn, disallows neutralino masses < 15 GeV (see Fig. 2 



of Ref. |46[). At variance with the conclusions of Ref. 

. in the present paper it is shown that an appropri- 
ate MSSM scheme fits the DAMA and CoGeNT annual 
modulation results quite well in force of the properties 
spelled out in Sect. Mil 

Ref. [6^ refers to our approach in the analysis of the 
CoGeNT data as being "somewhat unphysical", since , 
according to Ref. 63j, we would accept negative back- 



grounds. This is manifestly not the case, as can be easily 
understood by means of Eq. {T]), which defines the sta- 
tistical estimator we use in our analysis. We explicitely 
enforce the bound arising from the total rate, in order not 
to accept modulation amplitudes which would be incom- 
patible with the measured total rate. The last term in 
Eq. ([I]) does, in fact, penalize the \ 2 when the calculated 
rate becomes exceedingly large, stastitically incompatible 
with the measured total rate. We therefore do not accept 
negative backgrounds, contrary to the claim in Ref. 63 [. 
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